Simple models for granular force networks 
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A remarkable feature of static granular matter is the distribution of force along intricate networks. 
Even regular inter-particle contact networks produce wildly inhomogeneous force networks where 
certain "chains" of particles carry forces far larger than the mean. In this paper, we briefly review 
past theoretical approaches to understanding the geometry of force networks. We then investigate 
the structure of experimentally-obtained granular force networks using a simple algorithm to obtain 
corresponding graphs. We compare our observations with the results of geometric models, including 
random bond percolation, which show similar spatial distributions without enforcing vector force 
balance. Our findings suggest that some aspects of the mean geometry of granular force networks 
may be captured by these simple descriptions. 
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I. INTRODUCTION 

Dense granular materials are composed of many par- 
ticles interacting through multiple, persistent contacts. 
Without knowing the history of a particular granular as- 
sembly, the mobilization of the friction of each contact 
is indeterminate Even without friction, the simplest 
granular materials, such as packings of macroscopically 
uniform spheres, can still display highly non- uniform net- 
works of contacts — any difference in the size of parti- 
cles, including microscopic aberrations on the surface of 
"smooth" spheres, is sufficient to create a complicated 
contact network [2]. 

Experiments show that for any given particle within 
a granular material, the fluctuations in the magnitude 
of the force carried by that particle can be extremely 
large 0, 0] . This is the result of long-range force corre- 
lations that can be characterized as a force network in 
which some bonds of the contact network carry forces 
much larger than the mean force on all contacts. As in 
Fig. [TJ2, forces propagate along roughly linear "chains" 
of contacts. Because the force network depends upon the 
contact network, the propagation direction and distribu- 
tion of the chains depends strongly upon the history of 
the granular material d, 0] . 

The overall distribution of forces in a force network 
is roughly exponential, indicating the possibility for ex- 
tremely large forces. Since tangential forces are sup- 
ported only by friction, their distribution falls mono- 
tonically, while forces along the chains are peaked at a 
mean value before falling exponentially [Tj]. As grains 
are added, if the force network cannot support the gran- 
ular material, the material will rearrange until the new 
contact and force network supports the granular mate- 
rial. For this reason it has been suggested that granular 
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materials exhibit self-organized criticality, in which the 
force network behaves as an order parameter in response 
to displacements of the material 2]. 

For any configuration of grains (e.g.. Fig. [TJl), if the 
position and shape of each grain is known then a con- 
tact network can be generated describing the contacts 
between the grains. The equations describing a granu- 
lar material can be grouped into the known equations 
describing the behavior of the individual grains and the 
unknown equations describing the contacts between the 
grains as a function of material assembly and history. 
These unknown equations can be cast as undetermined 
relationship between inter-grain forces and the geometry 
of the contact network Q . 

Attempts to resolve the indeterminacy of the consti- 
tutive relations have focused upon statistical approaches 
similar to thermodynamics. Ideally, it might be possi- 
ble to determine aggregate properties of granular mate- 
rials that do not depend upon knowing the specifics of 
the grains to arbitrary precision. Unfortunately there 
are two problems that have made such an approach non- 
trivial. First, a granular material at rest is only in a 
quasi-cquilibrium. An arbitrarily small addition of en- 
ergy can cause irreversible (plastic) rearrangement of the 
material. This is in part due to the second problem, 
the lack of an analog to temperature to "thermalize" the 
system. The actual energy due to temperature is far too 
small to move grains and, in the static case, there is no 
other quantity allowing the material to explore phase- 
space for lower energy states. 

Despite these complications, there is a rich tradition of 
models to understand both the distribution of magnitude 
and the s pat ial correlation of forces within granular mat- 
ter ITollll| . The applicability of these models to actual 
granular materials has been limited in part by being re- 
stricted to regular lattices of grains. Recent numerical 
work [H, however, suggests that the spatial dis- 

tribution of force within irregular arrangements of grains 
might indeed be characterized through comparison with 
models with regular geometries. 

Typically the force network is determined using a 




FIG. 1: Analysis of the two-dimensional granular assembly of 
bi-disperse polymer disks depicted in 1. The contact network 
for the disks can be thought as a series of vertices centered 
at each disk connected by edges wherever two disks are in 
contact. 2 Because the disks are photoelastic, when they are 
illuminated between crossed polarizers bright fringes appear 
in disks under stress allowing force chains to be clearly seen. 
A load has been applied to the top of the pile causing the 
creation of additional contacts and force chains supporting 
the load. 3 We examine the gradient of the stress image as 
a rough measure of the number of fringes in each disk, which 
is proportional to the applied force. 4 Finally, we apply an 
algorithm to the gradient image to find a graph describing the 
force network. 



threshold — for instance, only bonds carrying forces 
larger than some value, often the mean force, are con- 
sidered to be part of the force network. In numerical 
simulations, critical-like transitions from a sparsely dis- 
tributed network of force chains to a percolating network 
have been observed as the threshold criteria for includ- 
ing bonds is varied. In one case, the magnitude of the 
inter-grain mobilization of friction was used and in 
another that magnitude of inter-grain forces [l^l. It is 
important to note that these thresholded force networks 
will not be rigid, since they lack bonds that, while car- 
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FIG. 2: Left In the two-dimensional g-Model, a fraction q of 
the total weight on a particle on a regular lattice is transferred 
to one lower neighbor, and the remaining 1 — g is transferred 
to the other. As originally formulated by Coppersmith, et 
al., the values of q can be drawn randomly from different 
distributions. In the "critical" case, the value of q is limited 
to either or 1. Right To make rough analogy to the q- 
Model for more complicated geometries, we use a geometry- 
based g-Model in which the fraction of a vertex's weight that 
is distributed along a downward pointing edge is proportional 
to the sine of the angle that edge forms with the vertex. 



rying only small amounts of force, provide rigidity to the 
network. 

There have been many attempts to model the behavior 
of force networks, most famously the "g-Model" of Cop- 
persmith (1996). In the two-dimensional g-Model, de- 
picted in Fig. [21 granular materials are represented as a 
lattice of sites with two upward neighbors and two down- 
ward neighbors. Each site has an intrinsic weight that is 
added to any weight it receives from the neighbors above 
it. The total weight on a site is distributed between the 
two downward neighbors with one receiving a fraction q 
of the total weight and the other receiving \ — q. The 
q-Model allows the determination of the distribution of 
weight by working downwards from the top of the lat- 
tice. This process can be shown to be diffusive-like in 
the continuum limit Q. 

Diffusive systems, however, do not allow for long-range 
propagating forces like those seen in force networks. Ex- 
periments on small systems have shown that even averag- 
ing over multiple force chain networks does not support a 
diffusive description of granular materials [l^- However, 
in the case where q is allowed to only take the values of 1 
or 0, the q-Model does show chain-like clusters of bonds. 
This case will be referred to as the "critical g-Model" [§] . 

In all cases the g-model is a scalar model in that the 
distribution of weight, not force, is being described. The 
forces exerted upon a particular lattice site do not need 
to balance, and this in part explains the lack of chain- 
like structures in all but the critical case. Vector models 
that attempt to include force balance also exist. One of 
the simplest lattice models that includes force balance is 
the "tripod" model that, in two-dimensions, assigns each 
site three downward neighbors, one directly below and 
two offset to either side by a fixed angle. In the con- 
tinuum limit this model can give wave-like equations for 
propagating stress [l^. Additionally, other models have 
included torque-balance on each site, giving meaning to 
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FIG. 3: In our experimental silo for observing granular force 
networks bi-disperse disks made of a photoelastic polymer are 
confined vertically to two-dimensions between clear plastic 
sheets and aluminum spacers serving as sidewalls. The grains 
rest upon a piston that is lowered slightly to make uniform 
the mobilization of friction at the walls. The entire silo is 
illuminated between crossed polarizers, which allows fringes 
induced by applied stress on each grain to be imaged with a 
CCD camera. 



the finite extent of the grains [ll| . Including both force 
and torque balance allows derivation in simple cases for 
equations describing force equilibrium as functions of the 
contact distribution 1 171. 



II. GEOMETRIC APPROACHES TO 
GRANULAR MATTER 

Research into the general behavior of force networks 
in granular matter has shown that certain properties, 
including the exponential distribution of large forces, 
should be expected to arise regardless of the grain-scale 
model. In this case, some properties of granular materi- 
als may be phenomena of the statistics of force networks, 
and not the micro-mechanics of the grains [l^, [l^ . 

This insight has led to approaches to granular mate- 
rials from the mesoscopic level. By assuming the pres- 
ence of force chains and describing their behavior, it may 
be possible to avoid some of the complexity of g ranular 
materials and still achieve meaningful results [201 l2l|. In 
particular, network approaches to granular materials may 
bear fruit because the statistical mechanics of networks 
in general are beginning to be understood [23 |. 

The critical g-model is similar to random bond perco- 
lation on a lattice. Percolation theory provides a frame- 
work for understanding how small clusters of connected 
sites aggregate into larger clusters as a function of the 
probability that a bond is selected. There are well known 
values, depending upon the geometry of the lattice, for 
transitions from loose aggregates of clusters to system- 
spanning, "infinite" clusters [23|. In the case of the crit- 
ical q-model, bonds are selected at random from a set 



of edges representing the connections between the lat- 
tice sites with the restriction that each lattice site has 
only one load-bearing downward bond. The set of load- 
bearing bonds can be thought of as a subset of edges that, 
with the vertices as the lattice sites, forms a sub-graph 
of the contact network. Since weight at a given vertex 
is equal to one plus the size of the tree rooted at that 
vertex, the distribution of cluster sizes gives bounds on 
the distribution of weights 

In disordered granular materials the network of inter- 
grain contacts may perhaps serve in place of a fixed lat- 
tice for percolation. Taking the contact network to be 
a graph, then the force network may be thought of as 
the sub-graph consisting of edges with greater than some 
threshold force. 

Properly selecting the edges that comprise the force 
network requires consideration of vector force balance for 
each particular network. However, recent work has sug- 
gested that granular force networks belong to a specific 
universality class that can be characterized through anal- 
ogy to random bond percolation, allowing the distribu- 
tions of forces in a wide range of isotropic and anisotropic 
gra nular s y stems to be scaled onto a single distribution 

111, [alii. 

In this paper we make experimental observations of 
the force networks within a granular system and use a 
simple approach to extract the statistical distributions 
of graph properties for these force networks. We charac- 
terize force networks both in terms of the graph-theoretic 
properties of degree distribution and cluster-size, as well 
as g-model inspired properties, including the general- 
ized q-model weight and the size of the cluster resting 
on a given grain. We then compare our experimentally- 
observed graphs with those generated by simple geomet- 
ric models to determine in what ways simple statistical 
approaches do and do not describe the magnitude and 
spatial distribution of forces within dense granular mat- 
ter. 



III. METHODOLOGY 

We use the experimental setup depicted in Fig. [3] to 
observe force networks within actual granular matter. 
We confine 650 bi-disperse (one quarter diameter 0.9 
cm and the rest 0.75 cm) photoelastic polymer disks 
within a vertical, two-dimensional silo. When illumi- 
nated between crossed polarizers, photoelastic materials 
show bright fringes of light in proportion to the stress 
applied to the material. By taking the gradient of an 
image of a photoelastic disk along multiple directions we 
effectively "count" the number of fringes and, using an 
empirical, linear calibration, can obtain the force distri- 
bution throughout the granular material (see Fig. [T]3). 

Using the gradient image of a particle force network 
we convolve with an image representing a single stressed 
grain to locate the positions of all stressed grains in the 
image. We then assign an edge between all vertices sep- 
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FIG. 4: Illustration of the graph quantities measured. The 
cluster size s = 10 includes all the vertices. The size of the 
upward tree rooted at the white circle (indicated by directed 
edges), t = 4: includes all vertices connected by upward edges 
only. The degree of the white node is fc = 3. The numbers in- 
dicate approximate geometric g-Model weights w, determined 
for each vertex by assigning each an inherent unit weight and 
distributing the weight among all downward edges in propor- 
tion to the sine of the angle each edge makes with the vertex. 



arated by a distance in the range of 1/4 to I where / is 
a manually-adjusted parameter. The lower range is nec- 
essary because some fringe patterns within a single disk 
produce multiple points when convolved, which has the 
effect of artificially increases clustering of the graph. 

Adding an edge to any grain that already has the 
threshold value k^ax edges causes the longest edge to 
that vertex to be removed. We use the smallest kmax 
that does not significantly alter P{k). 

The specific values of the network-finding parameters 
used depends upon many properties of the images of force 
networks, including the image resolution and dispersion 
of the disks. These parameters, in addition to the correla- 
tion threshold for locating vertices, must be tuned to pro- 
duce force network graphs that accurately describe force 
networks. We assigned edges within the range I = 13 
pixels and we allowed a maximum of kmax = 5 edges. As 
discussed below, the distributions reported here do not 
change unexpectedly with variation of correlation thresh- 
old. The set of vertices and edges found this way are the 
experimentally obtained force network (see Fig. [T]4). 

Using ensembles of force networks we can determine 
the probability distribution function for the average num- 
ber of contacts per vertex and number distribution of 
cluster sizes. We can also record the observed force at 
each vertex and create an ensemble probability distribu- 
tion for force. 

Figure S] illustrates several graph quantities that we 
measure. In addition to cluster size and vertex degree, 



FIG. 5: The distribution of vertex degree peaks at fc = 2. The 
lack of a sharp uptick at fc = 5 indicates that truncating at 5 
has not significantly affected the distribution. These results 
are from analyzing one hundred experimentally-obtained force 
networks under 56g load. 

we measure the size of the upward tree rooted at each ver- 
tex. The upward tree is a graph made by treating each 
edge as a directed edge oriented against gravity. The size 
of the upward tree corresponds to the number of vertices 
that are at least partially supported by a given vertex. As 
illustrated in Fig. [2l weights can be assigned to each ver- 
tex and propagated along each downward pointing edge 
in proportion to the cosine of the angle each edge makes 
with the vertex. This allows for a g-Model-like weight 
distribution to be generated for non-lattice geometries. 
For all of these quantities an aggregate probability den- 
sity function is calculated. 

IV. RESULTS 
A. Experimental Results 

We analyzed one hundred force networks obtained 
from experimental images of grains in a vertical silo sub- 
ject to a 56 g overload that created an approximately 
uniform distribution of force chains, as in Fig. [TJ2. Fig- 
ure[5]shows the distribution of vertex degrees and, as ex- 
pected for chain-like structures, vertices most often have 
two edges. From geometric considerations, the maximum 
number of contacts for a single disk in a bi-disperse mix- 
ture of disks with radii rsmaii and riarge is: 

7r/2/ arcsin [rsmaii/'^^l [r small ^larqe )] (1) 

In this case kmax — 6.85 and there is no probability 
of seeing a vertex with degree higher than six. Thus, 
as Fig. [5] shows, the imposed cutoff of kmax = 5 used 
in the edge finding routine does not significantly affect 
the distribution of P{k) since fewer than a tenth of all 
vertices have degree larger than 4. 
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FIG. 6: The distribution of Us, the number of clusters of size s 
divided by the total number of clusters, appears to be roughly 
linear over many decades of probability and two decades of 
cluster-size (corresponding to the system size) when binned 
logarithmically. A power-law distribution would support the 
idea of Roux et al. 12] that the force network may perhaps be 
considered a critical phenomenon. 



The distribution of the number of observed clusters 
of size s (as a fraction of the total number of clusters) 
appears to follow a power law, as indicated by Fig. [6l 
This behavior lends support to the suggestion of Roux, 
et al. that the force network may be an example of 
a self-organized critical phenomena. Fits were obtained 
for regressions to the data (neglecting the end points) for 
a power law (n^ = ois"^), the product of a power law 
and an exponential (rt^ ~ bis^'^e'^^'^), and a log-normal 
(ris = ^e'^^^^~'^^^ ) distribution. In all cases, the ra- 
tio of the variance predicted by the regression to the 
observed variance (i?^) is greater than 0.95, indicating 
likely fits [l^l . The fits for the power law-exponential and 
log-normal distributions, however, both produce positive 
curvature indicating increasing likelihood with very large 
cluster size. We believe such distributions are unphysical 
since increasing probability with increasing cluster size is 
not normalizable, so that a power law with an exponent 
a2 = —1.6 seems most likely. 

As shown in Fig. [71 the number distribution of up- 
ward clusters is shallower (a2 = —1.0) than for whole 
clusters but still appears to be linear. We expect the 
distribution of upward trees to generally be shallower 
than the distribution of clusters since each large cluster 
contains many upward trees. Regressions to power law, 
power law-exponential, and log-normal distributions all 
have > 0.85, but only the fits for the power law- 
exponential distribution give the observed drop-off with 
large upward tree size. 

Figure [8] shows the probability distribution of inten- 
sity of the gradient squared image at the pixel location 
of each vertex. The intensity corresponds with the force 
on the disk at that vertex and thus Fig. [5] is an approxi- 



FIG. 7: The logarithmically-binned distribution of rit, the 
number of upward clusters of size t divided by the total num- 
ber of such clusters, appears linear on a log- log plot, similar 
to the cluster distribution in Fig.|S]but with a shallower slope. 
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FIG. 8: The force observed at each vertex computed from 
the gradient is non-monotonic and the greatest probability is 
for forces below the mean force in the network (indicated by 
F/{F) = 10"). 



mate force probability distribution. This non-monotonic 
distribution peaks below the mean force and then falls 
off faster than a power law. 

We then compare the actual force distribution to what 
would be predicted using our geometric g-Model, obtain- 
ing the distribution in Fig. [5] The (/-Model weights are 
sub-power law but nearly monotonic, unlike the actual 
force data. 

The geometric limit of 6 possible contacts per disk cor- 
responds to 6 potential edges per vertex, however the 
average number of contacts (coordination number) for 
a disordered granular material is 4 0]. Assuming that 
there are four contacts per vertex, we can calculate a 
chance p that a bond carrying a large force is present. 



FIG. 9: The probability density function of weight assigned 
to vertices in the experimentally-obtained force networks by a 
method similar to the q-Model. Unlike the actual force data, 
the curve is monotonic except for a slight dip. 



For a particular vertex i, pi = ki/A and the average value 

P = = J2f=i A:i/4/650, where N is the total 

number of vertices in the force network graph, and 650 is 
the total number of particles. For the observed networks, 
p = 0.2214. 

Additionally, since we expect 4 contacts on average, 
Fig. [5] indicates that grains in the force network are typ- 
ically in contact with one or two grains not in the force 
network. 

For each vertex we can also compute the clustering co- 
efficient — how many edges are shared by the neighboring 
vertices (2) possible edges 25 1. Unlike a random network 
where any vertex may be contact with an unlimited num- 
ber of other vertices at any position, the choice of our 
radii and the quasi-two-dimensional geometry of our silo 
restricts the number of neighbor-neighbor contacts be- 
tween n neighbors to n — 1. We calculate c — 0.1219 for 
the observed force network. This amount of clustering 
is similar to that seen in scale-free networks that exhibit 
power law behavior, such as the World-Wide Web and 
power grids (26j . although since we have restricted the 
number of neighbor-neighbor contacts this is not a per- 
fect analogy. 



Finally, we examine the robustness of our results. To 
obtain force network graphs from our experimental data 
we convolved our images with a single fringe image ker- 
nel that does not necessarily match all arrangements of 
fringes. We used an acceptability threshold of 0.5, on a 
scale where 1.0 is a perfect match and 0.0 is the oppo- 
site of the kernel. As indicated by Fig. [TOl varying the 
acceptability threshold between 0.4 and 0.6 only slightly 
affects properties like vertex degree distribution without 
significantly changing the cluster size power law. 
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FIG. 10: Top Distribution of vertex degree for experimen- 
tally obtained force networks as a function of acceptability 
threshold for convolution with our fringe kernel. The dashed 
line indicates the distribution for the acceptability threshold 
of 0.5 used for the analysis in this paper. Lowering the thresh- 
old from 0.6 (thickest line) to 0.4 (thinnest line) in increments 
of 0.05 increases the number of vertices "found" and skews the 
distribution away from chain-like structures towards higher 
vertex degree. Bottom Despite the change in vertex degree 
distribution, the distribution of cluster size still indicates a 
power law with the same slope for all thresholds considered. 
The dashed line again indicates 0.5. 



B. Comparison with Spatial Models 

Having characterized actual force networks, it is possi- 
ble to make comparisons to simple models. Using an en- 
semble of realizations of the critical 5-Model on a 10 x 40 
vertex triangular grid (chosen to reflect the aspect ratio 
and roughly the number of vertices in our experiments) 
as in Fig. 1111 1 , it was possible to determine distribution 
functions for vertex degree, cluster size, upward tree size 
and weights assigned by the g-Model. 

In addition to the g-Model results, graphs were also 
constructed via bond percolation. Two regular trian- 
gular lattices were considered, one where in addition 
to downward- angled edges there were horizontal edges 
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FIG. 11: From left to right: 1 An instance of the q-Model in 
the critical case. 2 A random network (p = 0.25) composed of 
edges from a triangular lattice oriented horizontally. 3 A ran- 
dom network (p — 0.25) composed of edges from a triangular 
lattice oriented vertically. 



(see Fig. [TT]2), and one where there were vertical edges 
(Fig. [Tlj3). For each reahzation a subset of these edges 
were chosen with a probabihty p to create a randomized 
graph. With six possible edges per vertex, the triangular 
lattice is similar to the bi-disperse system though angular 
distributions are severely restricted. 

A range of different p values were examined for both 
triangular orientations. In addition to the p — 0.2214 
that was determined experimentally, we examined values 
of p = 0.15 and p = 0.3 above and below the estimated 
p. Finally, we examined p = 0.5, a value greater than the 
percolation threshold, Pc = 0.3473, for a system-spanning 
cluster in an infinite triangular lattice. 

We generated 5000 realizations each for the q-model 
and the two triangular lattice orientations. For the bond 
percolation graphs both vertex degree and cluster size are 
independent of whether the basis is oriented horizontally 
or vertically and so the joint ensemble of 10000 realiza- 
tions was analyzed. 

The distributions of vertex degree for the critical q- 
Model and triangular lattices with different bond prob- 
abilities are shown in Fig. [T^ Since the g-Modcl has 
to give chain like structures, it is not surprising to see 
that P{k) peaks at fc = 2, indicating that most elements 
are members of a single chain. For a random graph con- 
structed with the bond probability p — 0.2214 calculated 
from our experimental data, the distribution does not 
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FIG. 12: The probability density function for degree vertex 
for the critical g-Model and random triangular graphs. The 
critical g-Model is indicated by a dashed line. The bond per- 
colation results for the experimentally obtained p — 0.2214 
are indicated by □. Other bond probabilities, p = 0.15, 0.3 
and 0.5, are respectively indicated by lines of increasing thick- 
ness. While the g-Model by design gives the expected peak 
at fc = 2 — indicating chain-like structures — since in the 
critical case of the g-Model fc < 3, the best match to Fig. [S] 
appears to be p = 3. 



peak appropriately at A: = 2, but instead at A: = 1 in- 
dicating many pairs of bonds but no clusters. Instead, 
the peak of the distribution for p = 0.3 is most like the 
experimentally-observed distribution. The distribution 
for p = 0.3 over-represents single bonds when compared 
with experiments, but the distribution of larger k values 
matches well. 

Also plotted in Fig. [12] is the distribution for p — 0.15, 
below the measured p and not agreeing with the exper- 
imentally measured distributions. Finally, the distribu- 
tion for p — 0.5, shows too much skew toward large clus- 
ters. 

As indicated in Fig. [T31 the (/-Model gives a distribu- 
tion (with logarithmic binning) of cluster size that falls 
linearly (02 = —1.53) on a log-logplot until broadly peak- 
ing at roughly a quarter of the system size (100 vertices). 
Because we have chosen the critical case of the q-model 
to ensure long chains of contacts, it is not surprising that 
we see more large clusters. In contrast, the cluster size 
distributions for graphs constructed from triangular lat- 
tices shown in Fig.[T3]are monotonic and roughly linearly, 
with decreasing slope as p is increased until the critical 
threshold is crossed, after which the slope for small clus- 
ter sizes increases and a sharp peak at the size of the 
system appears, indicating infinite clusters (similar "bi- 
modal" distributions were observed previously bond per- 
colation simulations of much larger (10^) systems [27|). 
The slope for p = 0.3 is found to be most similar to that 
of the experimentally observed networks: 02 = —1.77 
with « 0.97 compared to 02 = —1.6 for experiment. 
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FIG. 13: The distribution of number of clusters of size s for 
the critical g-Model. On a log-log plot the distribution of 
smaller clusters is roughly linear (dashed line indicates slope 
— 1.53 found by fitting to linear region), but there is a peak 
at large cluster sizes since the critical g-Model requires that 
each site has a downward neighbor. 



Figure [T3] shows the distribution of "weight" assigned 
via the g- model for the critical g- Model graphs. The 
model distribution is more skewed toward large weights 
than the g-Model-like weight assigned to the experimen- 
tally obtained graphs. This is unsurprising since, by con- 
struction, the critical g-Model is a connected graph. The 
distribution of upward trees in Fig. (THlis roughly linear, 
with nearly the same slope (02 = —1.52) as cluster size. 

In Fig.[T7l we have the distribution of weights assigned 
by the geometric g-Model. None of the models showed 
the experimentally observed peak for actual force found 
in Fig. [5] Graphs composed from edges in vertically- 
aligned lattices — where an additional downward edge 
is possible — showed slightly broader distributions than 
those that were composed from horizontally-aligned tri- 
angular lattices. In the vertically-aligned lattices, the 
vertices with weight below the mean are distributed more 
broadly than in the horizontal case. For both alignments, 
the distribution for p = 0.5 was the most similar to the 
distribution for graphs obtained from experimental data 
— this values is much larger than the p measured exper- 
imentally. None of these distributions can be said to be 
a good match for the force or g-Model distributions for 
experimentally obtained graphs. 

As Fig. [18] shows, the distributions of upward trees 
for graphs below the percolation threshold are simi- 
lar. For p — 0.3 and below, the distributions are well- 
described by the product of a power law and an exponen- 
tial {R^ > 0.98) however the coefficients are very different 
from the experimentally-observed distribution of upward 
trees. For p = 0.5 the distribution is qualitatively similar 
to the experimental distribution in Fig. [71 but is not well- 
described by any of the distributions considered (power 
law, log- normal and power law-exponential). Generally, 




10° 10^ 10^ 10' 

Cluster Size s 



FIG. 14: The distribution of the number of clusters of size s 
for random triangular graphs. The bond percolation results 
for the experimentally obtained p = 0.2214 are indicated by 
□ . Other bond probabilities, p = 0.15, 0.3 and 0.5, are re- 
spectively indicated by lines of increasing thickness. For all 
cases there is roughly linear behavior for small to moderate 
(~ 100) cluster sizes. The transition value to universal clus- 
ters is Pc = 0.3473 for an infinite triangular lattice and as 
expected the slope varies non-monotonically, decreasing as p 
is increased until abruptly dropping at p = 0.5 where a second 
peak at large cluster size indicates that infinite clusters are 
forming. 

when there are vertical edges the distribution of cluster 
sizes is broader. This is reasonable since the addition of 
vertical edges makes large trees more likely. 

V. CONCLUSION 

We have characterized the graphs generated from ana- 
lyzing experimentally-obtained force networks in granu- 
lar material. Assuming that there are an average of four 
contacts per grain, we calculate that the probability that 
a contact carries a large force is p = 0.2214. We observe 
evidence for a power law scaling of clusters of bonds with 
large force, possibly indicating a simply-characterized un- 
derlying geometry. 

When we compare our experimental networks to bond 
percolation using random sets of edges from triangular 
lattices, we find that geometric properties of force net- 
works can indeed be matched depending upon the proba- 
bility of the presence of a bond p. We find that we cannot 
match all properties of the experimental networks using 
any one bond probability, including p — 0.2214. How- 
ever, the value of p = 0.3 worked well for both the vertex 
degree and cluster size distribution. 

A slightly larger value of p may be an artifact of under- 
sampling the number of minimally stressed particles that 
are hard to detect visually but do contribute stressed 
bonds. It is possible that further refinement of the graph 
finding algorithm is needed to get the correct value of p 
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q-Model Weight w/<w> 



FIG. 15: Probability distribution of weight assigned by the 
g-Model for the critical g-Model graphs. Because there are far 
more vertices with weights below the mean, the distribution 
of g-Model weights is different from both the actual force and 
geometric g-Model weight distributions for experimentally ob- 
tained force networks. 















o 














o ^ 




o 



10° 10^ 10^ 10^ 

Upward Tree Size t 



FIG. 16: Number distribution of upward tree size for the 
critical q-Model. 

for actual force networks. 

Because of obvious similarities, we have made compar- 
isons with the g-Model and attempted to assign weight 
to, and determine "force" distribution of, vertices in 
both g-Model and bond percolation graphs. None of 
the methods studied provide good agreement with the 
experimentally-obtained force distribution. This indi- 
cates that though graph theoretic approaches may in fact 
capture the geometry of force networks, the actual force 



distributions require more elaborate, vector analysis. 




"q-Model" Weight w/<w> 



FIG. 17: Probability distribution of geometric g-Model weight 
for graphs assembled from horizontally- (a) and vertically- 
(b) aligned triangular lattices. The results for the experi- 
mentally obtained p — 0.2214 are indicated by □. Other 
bond probabilities, p = 0.15, 0.3 and 0.5, are respectively 
indicated by increasing line thickness. The two alignments 
differ in the distribution of lower weights. In both cases, only 
the largest, p = 0.5, approaches the distribution observed for 
experimentally-obtained networks (Fig. |9]). 




Uoward Tree Size t 



FIG. 18: Distribution of number of upward trees of size t for 
bond percolation graphs assembled from horizontally- (a) and 
vertically- (b) aligned triangular lattices. The results for the 
experimentally obtained p — 0.2214 are indicated by □. Other 
bond probabilities, p — 0.15, 0.3 and 0.5, are respectively 
indicated by increasing line thickness. Only when infinite 
clusters are present (p = 0.5) do the distributions change 
significantly, in both cases giving a distribution that is linear 
over the same ranges as observed experimentally (Fig. ^ . 
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